Discrete breathers in a forced-damped array of coupled pendula: Modeling, 

Computation and Experiment 
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In this work, we present a mechanical example of an experimental realization of a stability reversal 
between on-site and inter-site centered localized modes. A corresponding realization of a vanishing of 
the Peierls-Nabarro barrier allows for an experimentally observed enhanced mobility of the localized 
modes near the reversal point. These features are supported by detailed numerical computations 
of the stability and mobility of the discrete breathers in this system of forced and damped coupled 
pendula. Furthermore, additional exotic features of the relevant model, such as dark breathers are 
briefly discussed. 



Introduction. Discrete breathers (also referred to as in- 
trinsic localized modes or ILMs) have been the focal point 
of numerous studies over the past two decades (for a re- 
cent review of the relevant activity see e.g. (H). Such lo- 
calized modes of nonlinear lattices have been found to be 
fairly generic and constitute a new paradigm in nonlinear 
science whose applications extend to physics and chem- 
istry, as well as to biology. Their original proposal in the 
broadly applicable context of anharmonic lattices [1, Q 
and subsequent theoretical analysis (including the rigor- 
ous proof of their existence under appropriate conditions 
has led to numerous experimental realizations rang- 
ing from optical waveguides and photorefractive crystals 
to micromechanical cantilever arrays and Josephson junc- 
tions, as well as in Bose-Einstein condensates and layered 
antiferromagnets, among many others [U, 

Typically, these ILMs come into two principal varieties, 
namely the on-site centered modes 0] and the inter-site 
centered ones [1]. In the simplest ones among the rele- 
vant models supporting such structures, either the for- 
mer, or the latter modes are stable, and this property 
is maintained when relevant parameters, such as the lat- 
tice spacing, are varied However, more recently it has 
been illustrated that stability exchanges are theoretically 
possible between these two types of modes JJ, which re- 
sult either in a true or approximate Q vanishing of 
the, so-called, Peierls-Nabarro barrier, namely the en- 
ergy difference between these modes. This vanishing is, 
in turn, relevant in that it enables the localized modes 
to move along the lattice without facing such an energy 
barrier, thereby enhancing their mobility 

Our aim in the present work is to demonstrate an ex- 
perimental realization of a setting which enables such 
stability reversals between on-site and inter-site modes, 
namely an array of pendula under the action of driving 
and damping. A detailed analysis illustrates such stabil- 
ity exchanges to occur numerically, and the experimental 
findings validate this result. Furthermore, enhanced mo- 
bility of the localized modes is observed experimentally 
and computationally near the exchange point and addi- 
tional theoretically interesting features of the model are 



predicted (in a regime that is not currently accessible ex- 
perimentally). The structure of our presentation is as 
follows. In the next section, we explain the experimental 
setup, followed by the theoretical/numerical setup. We 
close the paper with a detailed discussion of the results, 
as well as some potential future directions. 

Experimental Setup. We investigate a driven and 
damped chain of pendula each comprised of a rod of mass 
m and length L ending in a weight of mass A'l. They are 
suspended along a line by a taut piano wire stretched 
across a metal frame and connected via springs, as in 
[To| . The system is driven sinusoidally along a horizontal 
direction with an adjustable amplitude A and a tunable 
frequency lu^. The resulting dynamics of the pendulum 
chain is recorded by an overhead camera, and individual 
frames can be used to extract pendulum angles. The 
governing equations of motion for a pendulum in the 
chain take into account the detailed effects of driving 
and damping: 
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A2en + FloI cos(wdi) cos(6'„) = (1) 



with / = Mi^ + ImL^, lj^ = ) [nig^ + MgL) and 

F — — while A2 represents the discrete Laplacian. 
Here, g is 9.80 N/kg, and (3 denotes the torsional spring 
constant. The nature of damping in the system warrants 
closer inspection. Two types of friction are present in the 
experiment - an on-site and an inter-site one. The on-site 
component (of prefactor 71) comes from air drag as well 
as contact friction at the suspension points. The inter- 
site component (of prefactor 72) has its origin in the en- 
ergy dissipation due to the twisting of the springs, and it 
should be related to angle differences between neighbor- 
ing springs. This draws interesting parallels to the lattice 
turbulence models of [ll| and recent models of granular 
crystals jl^ . Equation ([J) can be non-dimensionalized by 
introducing the following parameters: ojb = — , C" : 
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FIG. 2; (Color Online). Top: 2-parameter planes for (left) 
fixed A — 1.12 cm and (right) fixed M = 1.8 grams. Full 
circles in the left panel correspond to the solutions shown 
in Fig. [1] Bottom: moduli of Floquet multipliers (indicating 
stability, when they are all smaller than 1) for (left) 1-site and 
(right) 2-site breathers. In this case, F = 0.06 and M = 1.8 
grams are fixed, while uih is variable. An exchange of stability 
bifurcation is observed for cu^, £ (0.804, 0.805). In both cases, 
the bifurcation arising for higher u>h is a Naimark-Sacker one. 



FIG. 1: (Color Online). Snapshot illustrating the experimen- 
tal setup and the observation of an inter-site breather (top 
panel), and comparison of experimental and numerical pro- 
files of stable inter-site breathers (middle panels) and on-site 
breathers (bottom panels). In all cases, A — 1.12 cm and 
u!h = 0.87. Circles represent the numerical results whereas 
the full lines with error bars correspond to the experimental 
profiles. 



and a — j^. 

This model is similar to tliat of , but here the driv- 
ing term depends on the angular displacement at each 
lattice site and on the driving frequency. The defini- 
tion of C suggests two experimental ways of adjusting 
this parameter. For gross adjustment, one can change 
the springs between the pendula modifying the torsional 
stiffness and coupling /3. A finer (and more convenient) 
way of varying C is via / and luq. This can be done most 
easily by changing the end-mass M which affects / and 
loq (see the respective definitions) . Note that by changing 
lvq, ujb and F are also affected. 

The on-site damping constant, 71, was measured using 
just one isolated pendulum, tracking its amplitude decay 
over time (which yields 71//). Incidentally, 71 was mea- 
sured for a number of endmasses, M, and found to be 
approximately independent of endmass for M < lOg. 

The parameter 72 was estimated as follows. Starting 
with a group of three pendula connected in a line via 
springs, the two outer ones are held fixed vertically, and 



the amplitude decay of the center pendulum is tracked. 
The exponential decay constant thus obtained is then di- 
vided by the angular velocity w for this mode (in the lin- 
ear regime) , which can be measured independently. This 
number corresponds to an effective a, and it is larger 
than ai = ^i/{Iujq). 

Table U gives the measured values of the experimental 
parameters. The driver amplitude, A(= 1.12cm), was 
selected to yield a value F of around 0.06. 

TABLE I: Experimental constants 



P (Nm/rad) 


71 (g cm^/s) 


72 (g cm^/s) 


L (cm) 


m (g) 


A (cm) 


0.0165 


284 


70 


25.4 


13 


1.12 



Theoretical/ Numerical Setup. We consider the 
breather existence and stability in the (non-dimensional) 
equations of motion. Our computations start at the anti- 
continuous (AC) limit (C — 0), using the two (differ- 
ent amplitude) attractors of an uncoupled pendulum. 
Breathers close to the AC limit only exist if there are 
two coexisting attractors of the single driven and damped 
pendulum - both phase-locked to the external driving. 
The highest-amplitude attractor will correspond to the 
excited breather site(s) whereas the smallest-amplitude 
one will be chosen for the low-amplitude background. 
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FIG. 3: (Color Online). (Left panel) Numerical (for F = 0.15) 
and (right panel) experimental examples of breather mobility. 
In the left panel, M = 3.8 grams and Wb = 0.66, while the 
right panel example is for M — 1.8g and uib = 0.87. The 
vertically stacked experimental profiles are spaced one period 
apart. 



Then, the breather will be a collective oscillation at the 
driving frequency, uji,. For the analyzed parameter range, 
the attractors oscillate in anti-phase. As a consequence, 
the excited site of the breather oscillates in anti-phase 
with respect to the tails, in agreement with the exper- 
imental observations of [lO]. Upon creation at the AC 
limit, such solutions have been continued for arbitrary 
C ^ 0, through a path-following scheme in the real space, 
as proposed in this method allows us to indentify 
both stable and unstable solutions. 

The linear stability of the ensuing discrete breathers 
is studied by means of a Floquet analysis. To this end, 
we add a small perturbation ^„ to the breather solution 
leading to the following equation for the perturbations: 

Cn + aCn + g{t; 0„(i))^„ - CA2^„ - aaAa^ - (2) 

with g{t;9n{t)) = cos6'„(t) — Fui^ cos{uJt,t) sin6'„(i) 

Stability analysis is performed by diagonalising the 
monodromy matrix A4, which relates the perturbation 
at i = to that at t = T: 

(l:;a^-(l:S 

The linear stability of breathers requires that the mon- 
odromy eigenvalues (Floquet m ultip liers) must be inside 
(or at) the unit circle; see e.g. [U dl for details. Due 
to the inter-site damping, only one property (implied by 
the real character of the monodromy) can be extracted 
for the multipliers, which is that they appear in complex 
conjugated pairs (A,A*). 

Results and Discussion. Let us start by comparing typ- 
ical experimental and numerical breather profiles. Sub- 
sequently, we will explore the regions in parameter space 
where the two types of breathers are stable, again com- 
paring theory and experiment. Finally, we briefly discuss 
additional interesting system features such as moving and 
dark breathers. 



In order to compare with experimental profiles, we 
choose a chain with 13 pendula and fixed-end bound- 
ary conditions. Fig. [T] details the experiment-theory 
comparison. The circles indicate computed values and 
the line shows the experimental profile; experimental un- 
certainties are indicated by the error bars. The uncer- 
tainty for angles higher than 100 degrees is ±3 whereas 
it is ±1 if the angles are smaller. An exchange of sta- 
bility is identified both numerically and experimentally: 
for small masses, 2-site breathers are stable whereas for 
larger masses, 1-site breathers become stable. 

The experimentally accessible parameters are M, Wb, 
and F (via A). To explore breather stability and symme- 
try more fully, we generate two-parameter diagrams for 
a larger-sized array. To this end, we have chosen a chain 
of iV = 41 pendula (so that boundary effects are of lesser 
importance). Fig. ^a.) illustrates the breather stability 
for fixed A = 1.12 cm (left panel where M and cjb are var- 
ied) . We observe that as M is increased, at a given Wb , an 
exchange of stability bifurcation is observed: the 1-site 
breather gains stability, and (practically) simultaneously 
the 2-site breather solution loses stability. Conversely, 
for a given Af , an exchange of stability can be induced 
by varying Wb. Above a critical value (« 0.875), the 
system experiences a Naimark-Sacker (NS) bifurcation 
(viz. generalized Hopf or Hopf bifurcation of periodic 
orbits) leading to the eventual breather disappearance. 
For driver frequencies cjb below a critical lower value, we 
observe the non-existence of any kind of breather (even 
unstable ones). Notice that the four dots in Fig. [H^a) 
correspond to the parameters of Fig. [T]) with M=Og and 
M=0.9g falling into the 2-site regime at ujh—Q-87, whereas 
M=1.8g and M=2.8g are in the 1-site regime. 

According to the figure, the stability exchange can also 
be initiated by varying the driver frequency at a fixed 
endmass, and this was also tested experimentally. For 
M=1.8g, the 1-site ILM is found to be stable down to 
Wb = 0.83. At that point, a transition to 2-site ILM ap- 
pears to be initiated. The resulting 2-site ILM is, how- 
ever, difficult to stabihze in the system of 13 pendula, 
and it usually decays within 5-10 periods. The observed 
transition frequency agrees reasonably well with the pre- 
dicted value of 0.805. At M=1.2g, the transition can be 
observed more clearly. The on-site breather prevails at 
Wb = 0.96. Below Wb = 0.91 the inter-site breather is 
stable; the predicted value is, in this case, 0.852. The 
NS-bifurcation at larger Wb is also clearly observed ex- 
perimentally but for endmasses above M=6.8g (see also 

m- 

In Fig. [DJb), a F-Wb plane for fixed M = 1.8 g is de- 
picted. Here we see that 2-site breathers are stable for 
small LJb with an exchange of stability bifurcation taking 
place in a narrow range of frequencies around lui, — 0.81. 
Below uJh — 0.65 breathers are NS-unstable. Notice, 
again, that for small values of F breathers can not be 
found. Figs. (H^c) and (d) show a typical parametric 
dependence of the Floquet exponents for two different 
configurations. 
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FIG. 4: (Color Online). (Top panels) Dark breather profile for 
C = 0.01 and ojb = 0.9 (left) and cJh-C plane for the values 
specified in the text. (Bottom panels) Example of pseudo- 
dark breathers with C — 0.01 which are stable in the range 
where "real" dark breathers are exponentially-unstable [(left) 
Wb = 0.8 and (right) ujb = 0.86]. 



The above stability exchanges between the two prin- 
cipal configurations (and the corresponding vanishing of 
the energy barrier between them) renders credible the 
scenario of breather mobility near such parameter sets 
[7, 8]. Such mobility has been observed (even in fairly 
large lattices with w 100 sites) and typically requires ap- 
propriately large values of F. A typical example of such 
a traveling breather for F — 0.15 is depicted in Fig. [3] 
(wb — 0.66 and M = 3.8 g) Such observations are in 
agreement with experimental indications of mobility (al- 
beit admittedly in the fairly short experimental lattice), 
as shown in the bottom panel of the Fig. [3] The experi- 
mental parameters are M=1.8g (C = 0.8) and ujt, — 0.87. 
In this regime, a stationary ILM can be made to move 
easily by bringing a fixed boundary into the wings of 



the ILM; in response, the ILM will move away from the 
boundary. The figure depicts an ILM repelled from the 
left boundary; it moves fairly quickly at first - alternat- 
ing between 1-site and 2-site symmetry - then slows down 
and finally comes to rest three sites over. To demonstrate 
another interesting feature of the model, let us now in- 
vestigate the existence and stability of dark breathers. 
In Ref. the existence of dark breathers in the dissi- 
pative Klein-Gordon lattice with non-parametric driving 
was demonstrated. Figured] shows the existence and sta- 
bility range of dark breathers. As shown, the range of 
stability of dark breathers is somewhat limited. We have 
varied parameters Wb and C for fixed 7 = 0.001, 72 = 
and F = 0.10. For small Wb, the dark breathers are 
exponentially unstable and lead to dark-localized struc- 
tures. For Wb above the range shown in the figure, the 
breathers experience a NS bifurcation. This instability 
is, however, small and system-size dependent. It was 
numerically found that dark breathers do not exist for 
Wb ^ 0.95, while for C > 0.012, stable dark breathers 
also do not exist - a constraint which in our present ex- 
perimental setting inhibits the formation of such coherent 
structures. 

Conclusions. We have illustrated that the setting of 
damped and driven coupled pendula is an ideal play- 
ground for showcasing many interesting features of dis- 
crete breathers. In particular, in addition to bearing un- 
usual model features, such as inter-site dissipation, this 
system exhibits complex stability properties evidenced 
by stability exchanges between on-site and inter-site con- 
figurations. These, in turn, lead to interesting breather 
dynamics, including the potential for mobility of such lo- 
calized excitations. Our numerical results based on the 
proposed model were found to be in excellent qualitative 
and good quantitative agreement with the experimental 
findings. Determining whether multi-site excitations or 
more exotic features such as dark breathers are possible 
would be among the natural extensions of the present 
work. 
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